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ABSTRACT 

Gas-giant planets that form via core accretion might have very different characteristics from 
those that form via disk-instability. Disk-instability objects are typically thought to have higher 
entropies, larger radii, and (generally) higher effective temperatures than core- accretion objects. 
In this paper, we provide a large set of models exploring the observational consequences of high- 
entropy (hot) and low-entropy (cold) initial conditions, in the hope that this will ultimately 
help to distinguish between different physical mechanisms of planet formation. However, the 
exact entropies and radii of newly- formed planets due to these two modes of formation cannot, 
at present, be precisely predicted. It is possible that the distribution of properties of core- 
accretion-formed planets and the distribution of properties of disk-instability-formed planets 
overlap. We, therefore, introduce a broad range of "Warm Start" gas-giant planet models. 
Between the hottest and the coldest models that we consider, differences in radii, temperatures, 
luminosities, and spectra persist for only a few million to a few tens of millions of years for 
planets that are a few times Jupiter's mass or less. For planets that are ~five times Jupiter's 
mass or more, significant differences between hottest-start and coldest-start models persist for 
on the order of 100 Myrs. We find that out of the standard infrared bands (J, i?, K, L' , M, 
N) the K and H bands are the most diagnostic of the initial conditions. A hottest-start model 
can be from '^4.5 magnitudes brighter (at Jupiter's mass) to ^^^9 magnitudes brighter (at ten 
times Jupiter's mass) than a coldest-start model in the first few million years. In more massive 
objects, these large differences in luminosity and spectrum persist for much longer than in less 
massive objects. Finally, we consider the influence of atmospheric conditions on spectra, and 
find that the presence or absence of clouds, and the metallicity of an atmosphere, can affect an 
object's apparent brightness in different bands by up to several magnitudes. 



Subject headings: planetary systems 
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1. INTRODUCTION 

Objects in the giant planet and brown dwarf mass 
ranges probably form in at least two main ways. 
On theoretical grounds, core accretion seems nec- 
essary to explain planets like Jupiter and Saturn 
([Pollack et al.l 1X9961) , while direct gravitational col- 
lapse via disk instability might be unavoidable for 
objects that form more t han several tens to hundreds 
of AU from their s tars (|Rafikovl [2007t iKratter et all 
[2OIOI : lRafikovl[20ll . If the different modes of forma- 
tion lead to different initial planet properties (such as 
heat content, radius, temperature, etc.) then there 
should be observables that, in principle, could be used 
to determine the n iode of formation of giant planets 
and brown dwarfs (iFortnev et al.l [20051 : iMarlev et al.l 
[2OO7I : iFortnev et al.ll2008D" 

Ambitious observing programs for the next 
decade hope to distinguish between core-accretion 
and disk-instability objects. Examples of such 
campaigns include the G e mini Planet Image r 
(GPI, iMacintosh et all [20061 : IMcBride et all [20ll . 
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the Near-Infrared Coronagraphic Imager (NICI) 
on Gemini South (|Liu et al.l l2010f ). the Spectro- 
Polarimetric High-contrast Exoplanet REsearch in- 
strument (SP HERE) on the Very Large Telescope 
([Beuzit et al.i 2008) , the Strategic Exploration of Ex- 
oplanets and Disks (SEEDS ) experiment on Subaru 
(|Tamural2009l:lJanson et al. 2011b ), and Project 1640 
on Palomar ([Hinklev et"ani2008L I2011[ ). When low- 
mass stars and brown dwarfs occur in close bina- 
ries, dynamical mass estimates can be derived, which 
can then be used to inform the initial conditions 
empl oyed in theoretical models (jKonopackv et all 
[2OIOI ). It might be very difficult to obtain similar 
mass estimates for directly imaged (wide-separation) 
planetary-mass binary companions to stars. Never- 
theless, spectra or multiband imaging of young plan- 
ets, together with estimates of the ages of the systems 
in which they are found, might be sufficient to derive 
constraints on both masses and initial conditions si- 
multaneously. 

Several authors have explored ways to observa- 
tionally discriminate b etween different formation sce- 
narios. IMarlev et all (2007, hereafter M07) exam- 
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ined the photometric signatures of "hot start" and 
"cold start" formation scenarios. In their framework, 
the hot-start objects represent disk-instabihty plan- 
ets and the cold-star t ones represen t thos e that form 
via core accretion. iFortney et al.l (2008) examined 
the spectral signatures that distinguish this partic- 
ular set of "hot-start" and "cold-start" models, and 
the influence of metallicity on emergent spectra.^ 

In this paper, we build on previous work by ex- 
amining the photometric and spectral signatures of 
a wide range of initial conditions, and the influence 
on spectra of different atmosphere types. Because it 
is not clear what precise initial conditions ought to 
be expected for either the core accretion or the disk 
instability mode of formation, we introduce a class of 
"warm-start" models. It is possible that both the core 
accretion and the disk instability modes of formation 
produce initial conditions that span a range, with disk 
instability forming planets generally - but not always 
- larger, hotter, and having higher entropies. Our 
warm-start models, then, might be representative of 
possible outcomes of either mode of formation. In ij2l 
we describe the dependence of planetary radius on 
mass and entropy, and the range of expected initial 
radii and entropies as a function of mass. In ij3l we 
describe the cooling and shrinkage of young objects 
with a range of masses and initial entropies. In ^ 
we discuss the evolution of hot- and cold-start planet 
models. In ij5l we present a continuum of warm-start 
models and discuss how observations could be used 
to constrain mass, initial entropy /radius, and atmo- 
sphere. Finally, in [JBl we summarize our conclusions. 



2. FORMATION SCENARIOS 

There are two main ideas about how objects form 
in the mass range of giant planets and brown dwarfs 
(from a fraction of Jupiter's mass - Mj - to tens of 
Mj)."^ One is that solid cores, similar to, or some- 
what more massive than, the terrestrial planets, form 
first, and then a runaway process of gas accretion fol- 
lows. In this picture, once a solid core of ■^lO Earth 
masses has developed, the rate at which it accretes 
"atmosphere" increases dramatically. Surrounding 

^ An alternative potential observational discriminant be- 
tween formation scenarios is related to the metallicity of stars 
where directly imaged giant planets are found. Giant plan- 
ets in the inner few AU around their stars presumably formed 
via core accretion, and these objects seem to o ccur more frc- 
quently around mor e metal-rich stars dFischer fc Valenti,2005i : 
I Johnson et al.ir2010l) . ICrepp fc JohnsonI l l201ll) . therefore, sug- 
gest that finding a similar correlation between planet occu- 
rance and metallicity for more distant (directly imaged) ob- 
jects might provide evidence that they form via core accretion, 
as well. 

^ Planets are often distinguished from brown dwarfs via 
the deuterium-burning criterion, often sai d to be --^IS Jupiter 
masses (Burrows et al. 1997, 2001; Chab rier fc Baraffd l2000l : 
[Chabricr et al. 2005). However, this mass depends on metal- 
licity, helium fraction, the presence or absence of clouds etc. 
(Saumon & Marlcvi i200S : SDicgcl et al. 2011). Mode of for- 
mation is arguably a better, albeit less readily o bservable, 
taxonomical discriminan t |[Bumjws_et_aLi2001i; .Chabrier et al.l 
l2005l : ISpiegel et aIll20Tll : [Schneider eTaTll25Tl] ) . 



gas from the disk flows inward and forms an accre- 
tion shock near the boundary of the growing planet. 
This shock can be quite luminous (jHubickvi et al.l 
120051: iMarlev et aIll2007D . and might allow the in- 
falling gas to radiate away much of it s initial heat. 
Some recent observational work (e.g., iJanson et al.l 
l2011al ) suggests that core accretion is indeed the pre- 
dominant mode of planet formation. The other idea 
is that the protoplanetary disk becomes gravitation- 
ally unstable and collapses into fragments, forming 
planets directly. The former idea, called "core accre- 
tion" or "nucleation instability " was s uggested more 
than three decades ago (.Harria 119781 : iMizuno et al.l 
119781 : iMizunoiriQSOD . and has been further dev eloped 
more recently (e.g., Ilkoma et al] 120001 I2001D . The 
second idea, called " disk instabilit y," has existed for 
at least six decades (iKuipedFlQSH ). In recent years, 
this mechanism, which remains controversial in the 
inner ~10 AU of a planetary syste m (Rafikov 2005), 
has been advocated largely bylBossI f|2000>,2007M201Q) 
and lMaver et~aI1 (12001 l200l . 

One obvious difference between the mechanisms 
is that core accretion involves accretion onto a 
core. This could lead objects that form through 
core accretion to have higher bulk metallicity than 
those that form through disk instability. Higher 
metallicity has two competing effects. For a 
given bulk entropy, it causes an object to have 
a smaller radius ( Burrows et all [2000t iGuillot et al.l 
120061 : iLeconte et al, .2009; .Ibgui et al.l l2010( ) How^ 
ever, it also leads to greater atmospheric opacity, 
which retards the rate at which an object cools and 
loses its initial entropy (Burrows et al. 2007). 

A second difference between core accretion and 
disk instability has to do with the timescale of for- 
mation. Core accretion involves two distinct steps - 
the formation of a rocky core, and the subsequent 
accretion of gas. Simulations suggest that these pro- 
cesses tak e on the order of se v eral to tens of million s 
of years ()Pollack et al.l 119961: iHubickvi et"all l2005f) . 
Core accretion simulations of the formation of Jup iter 
bvl Pollack eraLl (|1996f ) and lLissauer et al.l (|2009l) in- 
dicate that the formation timescale depends sensi- 
tively on the areal surf ace density (a j nit) of the pro- 
toplanetary disk. The iPollack et al.l ()1996D simula- 
tions indicate that the time required to form Jupiter 
is of order several million years (although for differ- 
ent values of (Tinit the formation timescale can vary 
by an order of magnitude). The opacities of grains 
in the protoplanetary nebula can also have an impor- 
tant influence on the core-ac cretion formation pro- 
cess (|Movshovitz et al.) I2010D . Direct collapse, on 
the other hand, is thought to proceed on an orbital 
timescale or a small multiple thereof - orders of mag- 
nitude faster than core accretion (Boss 2000). 

A flnal difference involves the heat of formation 
that remains in the object once it is fully formed. 
The speciflc entropy of gas in the protoplanetary disk 
is quite high compared with that of a bound ob- 
ject, such as a planet or a brown dwarf. But how 
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much of its initial entropy does the gas lose in be- 
coming part of such an object? In the core-accretion 
scenario, it is thought that an accretion disk forms 
with an accretion sho ck near the planet's boundary 
(jHubickyj et al.l [2005( 1. As gas flows through the lu- 
minous accretion shock, it lose s much of its initial en - 
tropy. For instance, in work bv lHubickvi et ahl ()2005D 
and M07, the gas is assumed to accrete at the entropy 
of the atmosphere of the forming protoplanet; i.e., it 
is taken to lose all its initial (excess) entropy. Core ac- 
cretion is sometimes, therefore, referred to as a "cold 
start" mechanism. In contrast, it is thought that disk 
instability leads gas to retain much more of its initial 
entropy, thus leading to the "hot start" appelation. 
For a given isolated object, a difference in initial en- 
tropy is observationally similar to a difference in age 
(the time required to cool from the high entropy con- 
dition to the low entropy condition) . If the age of an 
object is known precisely enough, therefore, it might 
be possible to distinguish a cold start from a hot start 
based on the apparent current entropy (and observ- 
ables that it influences, such as radius, effective tem- 
perature, luminosity, and spectrum). Since the rela- 
tive precision with which age must be determined to 
make such a discrimination increases with increasing 
age, younger objects are far more useful than older 
objects for testing formation models. 

The truth, however, might be somewhat more 
complicated than the simple "cold-start" / "hot-start" 
idea just presented. M07 pointed out that the ex- 
act post-accretion luminosity of core-accretion ob- 
jects depends on uncertain details of the accre- 
tion process; however, they emphasize that there 
should be a large separation between the initial en- 
tropies of core-accretion planets and disk-instability 
planets. This conclusion might turn out to be 
true, but is perhaps not yet definitively established. 
The accretion shock of the core- accretion scenario 
might no t be 100% efficient in rem oving th e initial 
entro py (jBromlev fc KenvonI 120111 : Morda sini et al.l 
120 111 ). Conversely, an accretion shock might form in 
the disk-instability scenario as well, radiating away a 
significant amount of an object's initial heat. In other 
words, the cold start might not be as cold as previ- 
ously thought, nor the hot start as hot. In principle, 
this could lead to a spectrum of "warm start" initial 
conditions. Nevertheless, there are generally thought 
to be some broad differences between the early prop- 
erties of objects that form via core accretion and those 
that form via disk instability. The infalling gas ini- 
tially has a significant excess specific entropy relative 
to the atmosphere of the forming planet onto which 
it is accreting. If this infalling gas retains some of its 
initial excess entropy, the initial entropy of a newly- 
formed planet would depend on this (unknown) frac- 
tion. If, on the other hand, the gas radiates away 
all its initial entropy down to the entropy of the gas 
that is already present in the protoplanet, the object's 
entropy at the conclusion of the accretion process de- 
pends on the (unknown) entropy of the gas that was 
already present. In either case, there is not an unam- 



biguous prediction for the initial entropy of a "cold 
start" object that forms via core accretion. 

In order to explore the range of initial entropies 
that are reasonable within the core-accretion context, 
we have computed a large number of new evolution- 
ary models. We emp loy the boundary conditions of 
iBurrows et al.l ()1997| ) , and we assume with M07 that 
the infalling gas radiates until it is at the same en- 
tropy as the gas already present. Figure [T] displays 
our new models and the MOT ones. The thick yellow 
bands in the top (entropy) and bottom (radius) panels 
represent so-called "hot-start" models, whose initial 
entropies range between those of M07 (open circles) 
and the mustard-colored curve (in the top panel) that 
represents our "hottest-start" model. The blue bands 
in each panel represent "cold-start" models, bounded 
by a pair of our new core-accretion models ("Acc: 
3x10-5" and "Acc: 3x10"^ Z: 10"^", respectively). 
The various new core-accretion models are described 
in 

2.1. The Influence of Radiating Efficiency During 
Formation 

Different model assumptions about the process 
of accretion infiuence the properties of newly-formed 
model planets. We, therefore, consider a range of 
assumed radiating e f ficienc ies of the accretion shock. 
iBromlev fc Ken-sroil ()2011t ) and Kenyon & Bromley 
(2011 in prep.) quantify this efficiency as l — r], where 
77 is the fraction of accretion energy that a parcel of 
gas retains when it hits the planet's photosphere, and 
Kenyon (2011, private communication) has provided 
us with several models of forming planets with var- 
ious values of 77. We estimate the specific entropies 
of these model planets by finding the values of spe- 
cific entropy that best explain (within the context of 
the Burrows et al. 1997 models) their combinations 
of mass, radius, and luminosity. These models used 77 
values of 0.15, 0.35, 0.50, and 0.75. The three models 
with a low value of 77 = 0.15 (which is, nevertheless, 
larger than the implied rj value of in both the M07- 
clone models and our newly calculated core- accretion 
models presented in Fig. [1]) span the mass range from 
0.43 Mj to 4.4 Mj, and for these, we infer specific 
entropy values from 5.2 to 6.3. These and all subse- 
quent specific (per-baryon) entropy values are in units 
of Boltzmann's constant per baryon. For compari- 
son, Jupiter has a specific entropy of ~6. Though it 
might initially be surprising that newly formed plan- 
ets should appear as cold as Jupiter is at ~4.5 billion 
years, the differenc es between the equation o f state 
(EOS) used in t he IBromlev fc KenvonI 1)20111 ) mod- 
els and the one (jSaumon et al.lll995D used in both 
the M07 models and in our own might expla in this 
apparent mystery. Since IBromlev fc KenvonI (|20ill ) 
assume that infalling gas retains 15% of its excess in- 
ternal energy, we presumably would have inferred a 
specific entropy closer to 9 (for reasons discuss ed be- 
low) if they had used the lSaumon et al.l (|1995D EOS. 

It is particularly instructive to examine the effect 
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of larger rj. Inferring the specific entropy from mass 
and radius as before, we estimate that a 0.88 Mj ob- 
ject formed with 77 = 0.35 has a specific entropy of 
7.5; a 1.2 Mj object formed with 77 = 0.50 has a 
specific entropy of 10.1; and a 1.6 AIj object formed 
with 7] = 0.75 has a specific entropy of 10.6. Given 
that the r] = 0.15 objects appeared to be roughly 
3 units lower in specific entropy than the cold-start 
band in Fig. [J (encompassi ng our own models and 
those of iMarlev et al.l l2007f ) , these specific entropy 
values at larger values of rj might similarly be reduced 
below what our models would calculate if embedded 
in the same formation scenario. Specific entropy is, 
of course, not linearly additive, and so the 10.6 value 
cannot be adjusted upwards in a simple way. Still, 
this suggests that if the gas cools fairly inefficiently 
while falling onto a forming giant planet, the final 
assembly might be significantly hotter than core ac- 
cretion models typically assume. 

2.2. Evidence for Hot Starts? 

Indeed, there is mounting evidence that suggests 
that particularly massive young objects form fairly 
hot. Several years ago, the brown dwarf eclipsing bi- 
nary 2MASS J05352184-0546085 was discovered. The 
two objects in this system have masses 57 Mj and 
36 Mj, and are thought to be quite young, less tha n 
a few million years old fStass un et all [2006I . [200I . 
Their radii (6.5i?j and 5.0Rj, respectively, where 
Rj is Jupiter's radius) and effective temperatures 
(2650 K and 2790 K, respectively) are indicative of 
the high-entropy initial conditions expected in the 
case of direct gravitational collapse. The masses and 
orbits of the objects in this system are more remi- 
niscent of very low-mass binary stars than of a star- 
planet system, and so the formation process might be 
expected to be more similar to star formation than to 
core accretion poergensll2608l ). 

Observations suggest that some lower mass ob- 
jects also form hot. For instanc e, the plan- 
ets found in the IIR8799 system ()Marois et al.l 
120081 120101 ) appear to be inconsistent with classical 
"cold start" models (Dodson-Robinson et alj I2009t 
Bowler et al. 2010; Janson et al. 2010; Cur rie et al.l 
20111: iMadhusudhan et al.l 12011 ). The planetary 
companion to /?-Pictoris (/3-Pic b) also appears to 
have formed ho tter than cold-start n iodels would 
have predicted (jLagrange et al.l 120101 : IQuanz et al.l 
[2OIOI) . Other young substellar objects, such as GQ 
Lupi b, can similarly inform our understanding of 
the character of their geneses and early e volution 
()McElwain et al.ll2007HNeuhauser et"aIll2008D . Given 
both the theoretical prejudice against disk instability 
as the dominant mode of giant planet formation and 
the recent observational evidence that objects do not 
form as cold as some core-accretion models suggest, 
the time is ripe to revisit what observations can tell 
us about the properties of very young planets. 

3. EVOLUTIONARY COOLING 



As a self-gravitating sphere ages, it cools and 
shrinks. The cooling is rapid at first, and slows as 
the object reaches lower effective temperatures. More 
massive objects retain their heat for longer, since their 
gravitational energy scales with (and by Virial 
equilibrium, the thermal energy scales with the grav- 
itational energy). 

Figure [5] vividly portrays how the cooling time 
scales with mass a nd with initial entropy, for 
iBurrows et al.l ()1997f ) models. For objects of 1, 3, 
10, and 15 times Jupiter's mass, we plot the time 
required to reach various threshold specific entropies 
(10, 11, and 12), as a function of the object's initial 
entropy. 15-Mj objects require several tens of Myr to 
cool to S* = 10 (where S is specific entropy, in units 
of Boltzmann constant per baryon), several Myr to 
cool to S = 11, and about 1 Myr to cool to 5* = 12. 
Less massive objects cool to 5 = 10 in a few Myr 
or less. Objects up to 10 Mj cool to an entropy of 
12 extremely quickly (in less than a million years af- 
ter "formation" - the end of accretion), irrespective 
of initial entropy. Reaching an entropy of 11 takes 
a factor of several longer (~2 Myr for 10-A/j models 
that start very hot). Reaching an entropy of 10 takes 
10- Mj models ~5 Myr. The cooling times are a fac- 
tor of ~4 less for 3-Mj models and another factor of 
^^3 less for 1-Mj models. 

As an object loses entropy, its radius shrinks. Fig- 
ure [3] demonstrates the dependence of radius on en- 
tropy for hydrogen-helium mixtures of mass 1, 3, 6, 

10, and 15 AIj. The radius increases rapidly with in- 
creasing entropy, and this rapid scaling is more dra- 
matic for lower mass objects. The intersections of the 
radius/entropy curves with vertical lines at S" = 10, 

11, and 12 show the radii that may be expected 
for (coreless) objects at the indicated entropies. At 
S — 12, radii of objects from 3Mj to 15Mj range 
from ~4 Rj to ~3 Rj. At S = 10, radii range 
from ~2 Rj (at IMj) to -1.5 Rj (at 15Mj). In an 
Appendix, we present an analytic fit to the radius- 
entropy-mass relation. 

4. "HOT START"/ "COLD START" SCENARIOS 

4.1. Initial Hot and Cold Start Conditions 

IMarlev et all ()2007D published a set of "Hot Start" 
and "Cold Start" models, where the initial entropy 
had a tuning-fork shape as a function of mass. In 
their framework, the cold-start models were intended 
to represent the result of a core-accretion process, and 
the hot-start models were intended to represent ob- 
jects formed via disk instability. 

We calculate a broad range of models, from ex- 
tremely hot (high-entropy, large radius) ones near 
the limit of gravitationally-bound stability through 
several relatively cold models that might represent 
possible outcomes of a core-accretion process. These 
models are summarized in Figs. [1] and |4l with the 
top panel of each figure showing initial entropy and 
the bottom panel showing initial radius as a function 
of mass. Figure [1] shows a total of seven families of 
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models. Two are the "hot" and "cold" start models 
from M07, shown with open and filled circles, respec- 
tively. Another is the "Hottest Start" family of mod- 
els, which represents an approximate upper limit to 
the initial entropy that might be expected of bound 
objects that form via a disk-instability scenario. The 
other four all represent core-accretion scenarios, with 
various initial entropies, accretion rates, and atmo- 
spheric properties. All these models assume that ac- 
cretion flows onto an initial protoplanet of low mass, 
with a given starting entropy. Accretion proceeds 
with newly-accreted material always at the specific 
entropy of the gas that has already accreted, which 
is itself slowly cooling in accordance with the atmo- 
spheric boundary condition. Our accretion models 
have the following properties: 

• The "Acc: 3 x 10~^" models begin at an entropy 
of 9.8, and accrete material at 3 x 10~^Mj/yr, 
simi lar to the assumed a ccretion rate in M07. 
The iBurrows et all ()1997[ ) atmospheric bound- 
ary condition is assumed. Forming a 1-Mj ob- 
ject takes '^3 x 10^ years. 

• The "Acc: 10"^" models are identical to the 
"Acc. 3 X 10~^" ones, except accretion proceeds 
at 1 X 10~^Mj /yi, a factor of 30 slower. 

• The "Acc: 10~^, lower S" models are identi- 
cal to the "Acc. 10~^" ones, except the initial 
entropy is 8.9 instead of 9.8. 

• The "Acc: 10"^, Z: 10"^" models are identical 
to the "Acc. 10^^" ones, except the atmosphere 
boundary condition is for 10^^-solar metallicity. 
The reduced atmospheric opacity allows cooling 
to proceed more rapidly. 

The entropy (top) and radius (bottom) panels of 
Fig. m also show several isochrones, corresponding to 
the "Hottest Start" models at various stages of evolu- 
tion, at 0, 1, and 10 Myr, and at 1 Gyr. Interestingly, 
the "Hott est Start" 1- Myr isochrone is quite simi- 
lar to the iMarlev et al.l hot start models ("M07 Hot 
Start"). 

Note that the IMarlev et al.l curves for hot and 
cold models nearly converge at IMj, and diverge pro- 
gressively more with increasing mass, leading to the 
tuning-fork shape seen in Fig. [TJ The lower track 
("MOT Cold Start") decreases with time because 10- 
Mj objects take longer to form than 1-Mj objects, 
so the already- ac creted g as has longer to cool. Al- 
though the Marle v~et al.l core-accretion models and 
the "Acc: 3 x 10~^" models nominally have the same 
accretion rate, the entropy of the Marlev et al.l mod- 
els clearly decreases with time, leading to the down- 
ward slope towards higher-mass objects (which take 
longer in forming), whereas the "Acc: 3 x 10"'°^" mod- 
els' initial entropy curve is nearly flat from 1-10 Mj; 
this difference might be attributable to the opacities 
and boundary conditions used. 



The upshot of Fig. [T] - particularly of the radius 
panel - is that there is a very dramatic difference 
in the initial properties of "Hot Start" and of "Cold 
Start" models. Newly-formed cold-start models are 
predicted to have radii in the ^ 1. 1-1. 7- i?j range, 
whereas newly-formed hot-start models might have 
initial radii of more than 3 Rj. Of course, since the 
core-accretion process might take of order millions 
of years and the disk-instability process is orders of 
magnitudes faster, the newly-formed objects under 
one scenario might correspond to a stellar age that 
is different by several million years from the newly- 
formed objects under the other forma tion scenario 
()Baraffe et al.|[2002l:lFbrtnev et al.|[2005l) . 

It is also worthwhile to keep in mind that planet 
formation models generally predict non-nionotonic 
evolution of luminosity during the formation process, 
with one or more large spikes in luminosity (e.g., 
IHubickvi et all [20051: iBromlev fc Kenvonll2011h . Our 
models begin (i.e., have t = 0) when the final mass 
assembly is complete. Although the very early evolu- 
tion of still-forming planets might contain interesting 
diagnostic observables, we do not explicitly calculate 
this process, and our models, therefore, do not pre- 
dict observables for this (comparatively brief) portion 
of planetary evolution. 

4.2. Evolution of Radius, Temperature, and Entropy 

We track the evolution of a broad range of mod- 
els, from "Hot Start" to "Cold Start," through an 
age of 1 Gyr. Figure U portrays our set of initial 
conditions (i.e., initial entropies and radii) for each 
mass. The large red and blue circles in Fig. |4] indi- 
cate our "Hot Start" and "Cold Start" models, which 
roughly trace the hottest and coldest starts shown in 
Fig. [TJ The medium-sized yellow and green dots rep- 
resent our variants of the M07 "Hot Start" and "Cold 
Start" models. 

Figure [5| depicts the evolution of specific entropy 
(top panel) , radius (middle panel) , and effective tem- 
perature (bottom panel) for our Hot Start and Cold 
Start models. Evolution curves are shown for a range 
of masses (1, 2, 5, 10 Mj). We also calculated the 
evolution of objects modeled after the M07 hot and 
cold starts (not shown). The qualitative behavior of 
the evolutionary trajectories of both our hot-/cold- 
start models and those of M07 is the same. At early 
times, objects that start hot are much larger and 
have much higher effective temperatures than objects 
of equivalent mass that start cold. As "Hot Start" 
(red) and "Cold Start" (blue) objects of the same 
mass evolve, their entropies, radii, and effective tem- 
peratures asymptotically converge. The convergence 
seen in Fig. [S] is more rapid at lower masses for two 
reasons: (1) lower mass objects evolve more quickly, 
and (2) our initial difference between hot and cold 
models is smaller for lower mass objects (similar to 
M07, though our hot / cold models encompass a some- 
what larger range). Within a few hundred million 
years, the memory of the initial conditions has been 
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lost even for 10- Mj objects. At early times, how- 
ever, the difference is quite stark for more massive 
objects (5, 10 Mj), where the initial radius and ef- 
fective temperature are each a factor of ~2 larger for 
"Hot Start" models than for "Cold Start" ones. These 
large differences in size (and, therefore, surface grav- 
ity) and effective temperature lead to striking differ- 
enc es in photometry an d spectra (as noted by M07 
and lFortnev et al.]|2008[ ). and are discussed in ij4.3l 

The difference between a "Hot Start" model and 
a "Cold Start" model may be thought of as a time 
difference (At) - the horizontal offset in Fig. [5] be- 
tween the red and blue curves of equal mass. As ob- 
jects reach ages comparable to or larger than this At, 
it becomes progressively more difficult to distinguish 
observationally between formation scenarios, on the 
assumption that these are the only formation condi- 
tions possible, an assumption that we revisit in [J5l 

4.3. Spectral Evolution 

The differences in entropy, radius, and effec- 
tive temperature between "Hot Start" and "Cold 
Start" models translate into differences in spectra and 
broad-band magnitudes. In order to compute spectra, 
we assume various atmospheres. In particular, we 
consid er four atmosphere types from iBurrows et al.l 
(|2011[ ): hybrid clouds at solar metallicity (our "fidu- 
cial" atmospheres); hybrid clouds at 3x solar metal- 
licity; cloud-free atmospheres at solar metallicity; 
and cloud- free at 3x solar metallicity. Our plan- 
ets are modeled as isolated objects, assumed to be 
in radiative equilibrium, and their emergent spectra 
are calculated wit h the line-by-line ra.diative transfer 
code COOLTLUSTY ((Hubenv et al.ll2003HBurrows et al] 
[2Q06I) . 

Figure [5] shows spectra as a function of mass and 
age for our fiducial (i.e., hybrid clouds at solar metal- 
licity) "Hot Start" and "Cold Start" scenarios. The 
red and blue (hot and cold) curves in Fig. [6] corre- 
spond to the large red and blue circles, respectively, 
in Fig. 21 Spectra are shown (assuming the source is 
at 10 pc) for objects of mass 1, 2, 5, and 10 AIj, and 
at ages of 1, 3, 10, 30, and 100 Myrs (spectra become 
progressively dimmer at later ages). At 1 Mj, hot- 
and cold-start models are nearly overlapping even at 
the earliest times. At higher masses, the strong dif- 
ference between the fiducial "Hot Start" and "Cold 
Start" scenarios becomes apparent. At ages of a few 
Myr, for instance, the 10-A/j objects are several or- 
ders of magnitudes brighter at short wavelengths and 
early times if they start hot than if they start cold. 
If the dichotomy between initial conditions is as dra- 
matic as in the fiducial models, then at early times 
there should be prominent observational differences 
(cf. M07: [Fortnev et al.if2Cl08l) . 

These differences in spectra imply differences in 
broad-band photometry. The evolution of J, H, K, 
L' , M, and iV-band absolute magnitudes (for hybrid- 
clouds, solar-metallicity models) is shown in Fig. [T] 
At young ages, the magnitude differences range from 



small ('^l) at IMj to more than 5 in some bands 
(e.g., J, H, K) at lOMj. As objects of a given mass 
evolve, their "hot-start" and "cold-start" brightnesses 
converge (in each band). The data in Fig. [7] are avail- 
able in Table [TJ and similar tables for the other three 
atmosphere types are available in electronic form on- 
line.'^ 

5. A CONTINUUM OF "WARM-START" PLANETS 

Planet formation theory informs our understand- 
ing of the inital conditions to be expected for core- 
accretion and for disk-instability scenarios. However, 
significant uncertainties remain. Consider again the 
range of hot-start and cold-start conditions repre- 
sented in both the entropy and radius panels in Fig.[T] 
A crucial assumption leading to the negative slope 
that is characteristic both of the IMarlev et al.l core- 
accretion models and of our own is that the accretion 
process adds no heat, meaning that each parcel ac- 
cretes only at the specific entropy of the gas that is 
already present. However, it is not clear that this 
assumption is warranted. Allowing infalling parcels 
to retain some of their nebular entropy would allow 
higher-mass core-accretion models to have higher ini- 
tial entropy than lower-mass models, thereby poten- 
tially filling in the "Warm Start" regions of the two 
panels of Fig. [1] 

Since we do not know a priori the initial entropies 
(and radii and temperatures) of objects that form via 
core accretion or via disk instability, it is useful to 
consider a broad range of possibilities. Indeed, if de- 
tails of the accretion process that might differ from 
one planet to another govern the fraction of nebu- 
lar entropy that is incorporated into the planet, then 
considering a broad range of possible initial entropies 
might be essential to reflect not just our current ig- 
norance, but also the range of actual initial condi- 
tions from which different newly-formed planets in 
our galaxy begin their evolution. 

Figure |4] summarizes our framework for treating 
a spectrum of initial conditions. At each of 1, 2, 5, 
and 10 Mj, we show a range of initial entropies (with 
different models spaced by 0.25 units): 8.75 to 10.5 
at IMj, 8.5 to 11.25 at 2Mj, 8.25 to 12.25 at 5Mj, 
and 8 to 13 at lOMj. As discussed previously (SHJ, 
the large red/blue circles indicate our hot /cold mod- 
els and the medium-sized yellow/green circles indi- 
cate the M07 hot and cold models. These particular 
choices for hot and cold models merely represent four 
values (at each mass) in a continuum of possibilities. 
The range of initial entropies in Fig. 2] corresponds 
to a range of initial radii: ^^1.4-3.3 Rj at IMj, ~1.3- 
3.0 Rj at 2Afj, ~1. 2-3.6 Rj at 5A/j, and ~1. 1-4.2 Rj 
at lOMj. 

5.1. A Dramatic Spread in Brightnesses and Spectra 

^ The model spectra presented in this paper are available in 
electronic form (at spectral resolution Ri204) at 

http:/ /www. sns.ias.edu/~dave/ 

and at[tittp://www.astro.prmceton.edu/~burrows/| 
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At a given age, models that start with different 
initial entropies can have very different brightnesses. 
For example, Fig. [S] shows the evolution of the K- 
band absolute magnitudes of the models denoted by 
dots in Fig.m At each of four different masses (1, 2, 
5, and 10 Mj), the evolution of if-band brightnesses 
is shown for a range of initial entropies. As seen pre- 
viously in Fig. [71 the memory of the initial conditions 
fades as objects evolve. At a given mass and age there 
is a range of possible K-hai\d absolute magnitudes, 
depending on the initial entropy. Cold-start and 
warm-start planets of higher-mass (>5 Mj) maintain 
relatively constant isT-band absolute magnitude dur- 
ing their early evolution. At 5 Mj, the near-constant 
brightness phase lasts for <10 Myr, and at 10 Mj, 
this phase lasts for several tens of Myr. Model plan- 
ets that begin with relatively high entropy are signif- 
icantly brighter in the early stages and evolve much 
faster, not experiencing the early-evolution plateaus 
in brightness that cooler-start models do. 

Figure [TU] summarizes the spread in _R"-band 
brightness as a function of initial entropy, initial ra- 
dius, mass, and age. In each panel, the ordinate indi- 
cates the increase in entropy relative to the coldest- 
start case (for the appropriate mass) depicted in 
Fig. 01 The corresponding initial radius scale for each 
is shown on the right y-axis as Ai?init , the increase in 
initial radius relative to the coldest-start case. The 
uneven spacing of Ai?init results from the fact that 
radius grows faster than linearly with entropy (see 
Fig. [3]). The color indicates the difference between 
the -ftT-band magnitude at the given entropy and the 
if-band magnitude for the coldest-start case. Cold- 
est start cases have initial specific entropies of 8.75, 
8.5, 8.25, and 8, and initial radii of 1.42, 1.32, 1.24, 
and 1.14 for masses of 1, 2, 5, and 10 Mj, re- 
spectively. At higher initial entropies, planets are 
from ~4.5 magnitudes (at Jupiter's mass) to ~9 mag- 
nitudes (at 10 Mj) brighter than the corresponding 
coldest-start cases. As planets age, the change in K- 
band absolute magnitude relative to the coldest start 
case approaches zero. More massive objects exhibit 
greater differences in brightness between hottest- and 
coldest-start cases, and maintain these differences for 
significantly longer. At 1 Mj, differences in bright- 
ness fade to less than a single magnitude by ^10 Myr, 
even for the hottest start. At 10 Mj, it takes several 
hundred Myr for the K-h&nd magnitude difference 
between a hottest-start and a coldest-start planet to 
reach 1 mag. Of course, at several hundred Myr, 
even 10- Mj planets are extremely dim in the K band. 
At longer wavelengths (e.g., M and N bands), cool 
objects such as ^1-Gyr old 10-Mj planets are still 
brighter than absolute magnitude 15 (not shown). 

It is instructive to examine the evolution of the full 
spectral energy distribution (SED) of model planets 
as a function of mass and initial entropy. Figure [Til 
shows the full spread in if-band absolute magnitude 
(from the hottest-start to the coldest-start models in 
Fig. [4]) for each of 6 bands, from J through TV. At all 
four planet masses shown (1, 2, 5, and 10 Afj), the 



difference in absolute magnitude is greatest in the K 
band (which is why we presented if-band evolution 
plots in Figs. [HI and [TU)) . The spread in magnitudes 
is smallest in M band, where at 1 Mj the difference 
between the hottest-start and coldest-start models is 
only «1.5 mags, even very early after "formation" 
(at 1 Myr) . The relatively small differences in bright- 
ness between hottest and coldest models at M and N 
bands are in large part because these bands are in the 
Rayleigh- Jeans tail of a planet's SED. As a result, the 
brightness in these bands is roughy proportional to ef- 
fective temperature, which varies by no more than a 
factor of a few at early times, and significantly less at 
late times (see Fig. [5]). 

5.2. The Effect of Atmosphere 

It is worthwhile to consider the effect of different 
atmosphere types (discussed in §4.3p on observables. 
Equal-mass objects with identical initial entropy can 
have very different colors depending on their atmo- 
spheric properties. Figure [51 provides an illustrative 
example of these differences, for a very hot start (ini- 
tial entropy of 13) lO-Afj model. The ordinate rep- 
resents the difference in magnitude, as a function 
of age, between each of four atmosphere types (hy- 
brid clouds, solar metallicity; hybrid clouds, 3x solar 
metallicity; cloud-free, solar metallicity; cloud-free, 
3 X solar metalhcity ; these are the four atmospheres of 
[Burrows et al.|[201ll ) and our fiducial atmosphere (hy- 
brid clouds, solar). In J band, differences start small, 
but by late times the cloud-free solar models are two 
magnitudes fainter than the fiducial models, and 3 
magnitudes fainter than the cloud-free 3 x solar mod- 
els. In K band, the differences between some models 
and others grow to more than 3 magnitudes. In A'' 
band, by contrast, the differences between cloudy and 
cloud-free models begin as moderate (~1 magnitude) 
and slowly decrease with age. 

Another perspective on the effect of atmosphere 
type is provided in Fig. I12[ which presents the evo- 
lution of the spread in absolute magnitude as a func- 
tion of atmosphere type and spectral band for 5 Mj 
planets. The four panels, representing the four dif- 
ferent atmosphere types, look broadly similar, with 
the largest differences between hottest and coldest 
planets generally occuring in H and K bands. How- 
ever, there are some quantitative and qualitative dif- 
ferences among the four panels. At solar metallicity, 
both the cloudy and the cloud-free models exhibit 
the largest hot-cold differences in K band ('^8 magni- 
tudes at early times), whereas at 3x solar metallicity 
H and K bands have approximately equal hot-cold 
differences (^7 magnitudes at early times). For all 
four atmospheres, M band shows the smallest hot- 
cold difference throughout the 1 Gyr of evolution dis- 
played. A key result shown in Fig. [T^l is that the 
effect of atmosphere type on the spread in absolute 
magnitude is not nearly so strong as the effects of 
changes in mass, age, or initial entropy, but spectra 
can nonetheless differ by several magnitudes depend- 
ing upon atmospheric properties. 
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One can view Figs. IWT^ in a forward sense. If 
we know the age, mass, initial entropy, and atmo- 
sphere of a model, these figures and associated tabu- 
lated data tell us the absolute magnitude in each of 6 
photometric bands in the infrared. However, this in- 
formation is probably more useful when viewed in an 
inverse sense. For comparing these models with data, 
the magnitudes or spectra will be the observed quan- 
tities, and one will try to make inferences about the 
properties of the object. These properties include the 
same parameters listed for the forward problem: age 
(about which there might be some prior constraints 
from either the star or its environment), mass, initial 
entropy, and atmosphere type. By checking which 
models reasonably match observed data, we might 
hope to find contours in a multidimensional param- 
eter space that rule in or rule out possible models. 
In this way, we may use data to infer the initial en- 
tropies of observed objects. Note that an exploration 
of our model space indicate s that the th i cker c louds 
employed in the models o f iCurrie et al.l ()2011[ ) and 
iMadhusudhan et all ()2011[ ) seem to be essential for 
finding self-consistent (coeval) models of planets b, c, 
and d in the HR8799 system. 

One contour that appears to be available for 
the H R8799 system is the stellar age. iMarois et al.l 
(|200^ argue that several lines of evidence sug- 
gest that t he system has an age in the range 30- 
160 Myr. IMova et al.l (j2010D find that asteroseis- 
mological data imply an age closer to 1 Gyr. How- 
ever, if the planets were 1 Gyr old, their luminosities 
would imply masses large enough to render the sys- 
tem dynamically unstable (iFabrvckv fc Murrav-Clavl 
20101: lMoro-Martmet"aII I2O1O0 . Furthermore, 
Zuckerman et al.l ()20H'i rfind that HR8799's proba- 
ble association with a nearby group suggests an age 
of ~30 Myr. We therefore assume, for now, that 
the ob jects have ages in t he range originally identi- 
fied bv lMarois et alj ()2008D . The fourth planet found 
in the system, planet e, was r ecently discovered i n 
Ki and L' ba nd observations bv lMarois et al.l ()2010D . 
IMarois et al.l find that the planet's luminosity sug- 
gests a mass of 7-10 Mj if its age is in the range 
30-60 Myr. We tried to fit the photometry (the ab- 
solute magnitudes in these bands are 12.93 ± 0.22 
and 11.61 ± 0. 12, respectively) with bot h the AE 
cloud models of IMadhusudhan et al.l (|2011[ ) and with 
our own set of models. Because of the paucity of 
data for planet e, a range of masses, ages, initial 
entropies, and atmosphere s are consistent with the 
available data. Both the Ma dhusudhan et al.l AE 
cloud models and our own hybrid cloud models (hyls 
and hy3s) are consistent with the data at th e ages 
and masses sugg ested by IMarois et ahl (|2010[ ). At 
these ages and masses, consistency with the photome- 
try requires significantly higher initial entropies than 
our coldest start scenarios. Since our cloud-free mod- 
els (cfls and cf3s) are dimmer in Kg and L' bands 
than the hybrid cloud models, fitting the data with 
the cloud-free models requires higher masses within 
the age-range (30-60 Myr) specified bv IMarois et al.l 



(j2010t ). Additional photometry at other bands could 
break some of the present degeneracies and might rule 
in or rule out some atmosphere types. 

6. CONCLUSION 

Gas-giant planets are thought to form in one of 
two scenarios - either runaway accretion onto a rocky 
or icy core, or direct gravitational collapse via a 
Jeans-like instability triggered in the protoplanetary 
disk. The former process is generally thought to lead 
to colder (i.e., lower entropy) objects than the lat- 
ter, but first-principles calculations cannot yet specify 
with certainty what the initial (post-formation) en- 
tropies of objects should be in the different formation 
scenarios, nor whether the range of entropies resulting 
from the different mechanisms might overlap. 

We seek observational determinants of the vari- 
ous models for giant planet origins, including their 
"initial conditions." Candidate properties that might 
be discriminating include metallicity, rotation rate, 
mass function, or orbits. Core-accretion-formed ob- 
jects might generally be expected to have higher bulk 
metallicity than disk-instability-formed objects; and 
the different modes of formation might tend to lead 
to differences in angular momentum, typical object 
mass, and initial orbit. A wrinkle associated with the 
idea of using rotation rate as a discriminant has to do 
with the fact that, for 10 Mj and slightly more mas- 
sive objects, a deuterium-burning phase occurs a few 
million years post-formation. This phase delays the 
evolutionary shrinkage, and when it has completed 
objects cool and shrink rapidly. It is conceivable that 
there is a signature of this rapid shrinkage embed- 
ded in the rotation rates of these objects, where mas- 
sive objects burning deuterium might appear to be 
rotating anomalously slowly, and conservation of an- 
gular momentum might cause a spin-up in the post- 
deuterium stages. But without knowing precisely the 
initial total angular momentum budget of a massive 
planet, it is unclear how such a signature could be 
interpreted. 

However, more direct signatures involve broad- 
band colors and spectra. We have presented a large 
suite of evolutionary models of non-irradiated gas- 
giant planets, which are publicly available in elec- 
tronic form. Our models span a range of "ini- 
tial conditions," with masses ranging from 1 to 15 
times Jupiter's, entropies ranging from 8 to 13 Boltz- 
mann constants per baryon, four different atmosphere 
types, and a dense range of ages up to a gigayear. 
Motivated by (i) the strong preference for the core- 
accretion mode of gas-giant planet formation within 
the inner tens of AU, (ii) the recent observations of 
several young exoplanets that seem to have been born 
hotter than traditional "cold start" models, and (Hi) 
the uncertainties inherent in the current generation 
of both core-accretion and disk-instability models of 
planet formation, we introduce a set of warm-start 
models that are at intermediate entropies between 
traditional "hot starts" and "cold starts." We find 
the following robust conclusions: 
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• Objects that begin at larger radius and higher 
entropy can be significantly brighter in the first 
few million years (by 2-9 magnitudes, depend- 
ing on mass and spectral band). 

• These photometric differences are most promi- 
nent in H and K bands, and least in M and N 
bands. 

• For lower-mass (1 and 2 Mj) gas giant planets, 
these differences fade to barely observable lev- 
els within ~10-30 Myr, whereas for more mas- 
sive (5-10 Mj) planets, a hottest-start model 
can remain several magnitudes brighter than a 
coldest-start model in the K band for more than 
100 Myr. 
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• Changes in atmosphere type can result in 
changes in brightness of a few magnitudes for 
some masses, ages, initial entropies, and spec- 
tral bands, but atmosphere type has less influ- 
ence on a planet's emergent flux than mass, age, 
and initial entropy. 



Some caveats remain. For one, our exploration 
of atmosphere types is far from complete, compris- 
ing only a single kind of cloudy atmosphere and a 
single cloud-free atmosphere, each at both solar and 
three-times solar metallicity. Other types of clouds, 
non-solar mixtures of elements, and nonequlibrium 
chemical processes could all have important affects 
on spectra that are not captured by our models. Per- 
haps more significant, we do not model in detail the 
formation processes, including the large spike in lumi- 
nosity that might occur at very early "pre-formation" 
times, due to an accretion shock during either direct 
collapse or core accretion. Furthermore, since disk 
instability is thought to form planets orders of mag- 
nitude faster than the '^several million year timescale 
for core accretion, there is a hidden age ambiguity 
that must be accounted for when our model planets 
are compared with observed systems. That is to say, 
perhaps a model that is supposed to represent disk 
instability should be taken to be several million years 
older than a model that is taken to represent core ac- 
cretion. Rather than try to build this into our public 
models, we prefer to simply start the models at time 
zero with a particular initial entropy and allow others 
to adjust the clocks as they consider appropriate. 

We hope that these models will prove to be useful 
in diagnosing observed young planets. By comparing 
data to model atmospheres, one can infer the initial, 
post-formation conditions of gas-giant planets. Even- 
tually, this type of analysis might constrain planet 
formation mechanisms by, for instance, quantifying 
how much entropy nebular gas retains as it accretes 
onto a core. In this way, models such as those in this 
paper contribute toward developing a coherent pic- 
ture of the evolution of planetary systems, from their 
formation through the epochs when we observe them. 
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Initial Entropy of Newly-Formed Planets 
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Fig. 1. — Models of specific entropy (fop) and radius (bottom) of newly-formed planets. Values from M07 are shown as open 
(filled) circles for hot (cold) start models. Four families of core accretion models are shown (accretion rates of 3 X 10~^Mj/yr to 
IQ-^Mj/yr, with various metallicities and initial entropies). There is also a "Hottest start" family of models shown with extremely 
high initial entropies and radii. 1-M yr, 10-Myr, and 1-Gyr isochrones of the "Hottest start" model are shown; the 1-Myr isochrone 
is fairly similar to the M arlev et al.l "Hot Start" models. Representative bands are shown for cold-start (blue), hot-start (yellow), 
and warm-start models. 
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Fig. 2. — Cooling time as a function of initial entropy for lBurrows et al] 1 119971 ) models, for 1, 3, 10, and 15 Mj. Line thickness 
indicates mass, with thicker lines denoting more massive objects. Red, magenta, and blue curves have final specific entropies of 12, 
11, and 10, respectively. Objects even as massive as 10 Mj cool to an entropy of 12 very quickly (in less than a million years after 
the end of accretion), regardless of their initial entropy. Cooling to an entropy of 11 takes a factor of several longer (~2 Myr for 
10- Mj models that start very hot). Cooling to an entropy of 10 requires ^^5 Myr for 10-A/j models. Cooling times are a factor of 
~4 less for 3-Mj models and another factor of ~3 less for 1-A/j models. 




Fig. 3. — Radius vs. Entropy for objects of mass 1, 3, 6, 10, and 15 Mj. Radius increases rapidly with increasing entropy, more 
so for lower mass objects. Vertical lines (in colors corresponding to those in Fig.|2)l mark where the entropy per baryon has values 
of 10, 11, and 12. 
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Fig. 4. — Range of initial entropies (top) and corresponding initial radii (bottom) that we employ for evolution calculations are 
displayed . For each of 1, 2, 5, and 10 Mj, a broad set of initial entropy values (in steps of 0.25) are listed. Very large red (blue) 
circles indicate our "Hot Start" ("Cold Start") scenarios. Medium-sized yellow (green) circles indicate our versions of the M07 
"hot-start" ("cold-start") models. For more massive objects, the highest and the lowest plausible entropies both become more 
extreme. The radii corresponding to the entropies of the top panel are shown in the bottom panel. 
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Fig. 5. — Evolution of specific entropy {top), radius {middle), and effective temperature {bottom), for a range of masses (1, 2, 5, 
10 Mj), for a particular set of "Hot Start" (red) and "Cold Start" (blue) initial conditions, corresponding to the large red and blue 
dots in Fig. |4l respectively. The differences in radius and temperature between the "Hot Start" and the "Cold Start" objects is 
dramatic at early times, particularly for more massive objects. As "Hot Start" and "Cold Start" objects of the same mass evolve, 
their entropies, radii, and effective temperatures asymptotically converge, such that within a few hundred million years the memory 
of the initial conditions has been lost even for 10-Afj objects. Lower mass objects lose memory of their initial conditions within a 
few to a few tens of million years. 
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Fig. 6. — Evolution of spectr a for "Hot Start" (red) and "Cold Start" (blue) models for objects with hybrid-cloudy atmospheres 
at solar metallicity (Burrows c t al.|[2011l) . The initial conditions correspond to the large red and blue dots in Fig. [3] Spectra are 
shown (assuming the source is at 10 pc) for objects of mass IMj {top left), 2Mj (top right), 5Mj {bottom left), and lOAf j {bottom 
right). In each panel, 5 spectra are shown for each initial condition, indicating isochrones of 1, 3, 10, 30, and 100 Myrs. At the top 
of each panel, thin black lines indicate the transmission functions for J, H, K, L' , M , and bands. 
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Fig. 7. — Evolution of absolute magnitudes. "Hot Start" (red) and "Cold Start" (blue) models are shown for 1, 2, 5, and lOAfj. 
Evolutionary trajectories are shown for J (top left), H (top right), K (middle left), L' (middle right), M (bottom left), and N 
(bottom right) bands. Hot and cold-start models begin at very different brightnesses, but as objects age the differences quickly 
approach zero. 
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FlG. 8. — Illustrative examples of evolution of the difference in magnitude relative to the fiducial atmosphere (i.e., relative to a 
hybrid clouds, solar metallcity atmosphere model). This figure shows the differences for an object of mass lOAf j with initial entropy 
of f3. In some bands, at some ages, the differences in brightness from one atmosphere model to another (at identical object mass 
and initial entropy) can be several magnitudes. 
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Fig. 9. — Evolution of absolute X-band magnitude for different masses and different initial entropies. Evolutionary tracks are 
shown for each of the models represented by dots in Fig. |4] (with corresponding colors). Thick red and blue curves represent our 
hottest- and coldest-start models, corresponding to the large dots of the same color in Fig.|4] medium-thick yellow and green curves 
similarly represent our version of the M07 "hot-start" and "cold-start" models. At each of 1 Mj {top left), 2 Mj {top right), 5 Mj 
{bottom left), and 10 Mj {bottom right), planets begin bright in the X-band and cool with time. Those that begin with higher 
initial entropy (and larger radius and higher effective temperature) are significantly brighter at early times than those that begin 
cooler. As objects evolve, the memory of the initial conditions fades, but at a given mass and a given age there are a range of 
possible i^-band magnitudes, depending on the initial entropy. 
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Fig. 10. — Evolution of if -band absolute magnitude as a function of initial entropy and mass. For each of four planet masses — 
IMj (top left), 2Mj {top right), 5Mj (bottom left), lOMj {bottom right) - the change in if-band absolute magnitude relative to 
the corresponding "coldest start" case in Fig. |4]is shown. The coldest start cases have initial specific entropies of 8.75, 8.5, 8.25, 
and 8, and initial radii of 1.42, 1.32, 1.24, and 1.14 Rj, for masses of 1, 2, 5, and 10 Mj, respectively. The left j/-axis of each panel 
shows the increase in specific entropy relative to the coldest start. The right j/-axis of each panel shows the change in radius relative 
to the coldest start. At higher initial entropies (and larger initial radii), planets are brighter in the K band. As planets age, the 
change in if-band absolute magnitude relative to the coldest start case approaches zero. Note that each panel has an independent 
color stretch. More massive objects exhibit greater differences in brightness between the coldest and the hottest starts, and their 
differences persist for longer. 
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Fig. 11. — Evolution of spread in absolute magnitude as a function of mass and spectral channel. For each of four planet masses 
- IMj (top left), 2Mj {top right), 5Mj (bottom left), lOA/j {bottom right) — the difference in absolute magnitude between the 
coldest start an d the hottest st a rt cas es in Fig. |4]is shown. The atmospheres of all objects are the hybrid cloudy models at solar 
metallicity from 'Burrows et al.l 11201 II ). Note that each panel has an independent color stretch. In all six bands, more massive 
objects exhibit greater differences in brightness between the coldest and the hottest starts, and their differences persist for longer. 
The spread in absolute magnitude is greatest at K band throughout the range of masses. 
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Fig. 12. — Evolution of spread in absolute magnitude as a function of atmosphere type and spectral channel. For each of four 
atmosphere types - cloudy at solar metallicity (top left), cloudy at three times solar metallicity (top right), cloud-free at solar 
metallicity {bottom left), cloud-free at three times solar metallicity {bottom right) - the difference in absolute mag nitude between 
the coldest start and the hottest start cases in Fig.[4]is shown. Details regarding the atmosphere types are described in lBurrows et all 
l(2011). All model planets are 5 Mj. The effect of atmosphere type on the spread in absolute magnitude is not dramatic, although 
significant differences exist from one spectrum to another which are not apparent in this figure. 
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APPENDIX 
AN ANALYTIC FIT FOR RADIUS 



Here, we present a fit to the radius-entropy-mass relation that works quite well for masses between 0.3 Mj 
and 15 Mj (and for reasonable radii). Because of a distinct change in the behavior of the R[S, M] surface near 
2 Mj, we present two separate fits in the two domain regions. Each is an exponential of a polynomial (3rd order 
in S and 2nd order in M). 

For M < 2 Mj, 

Poo = -1.27 
Pio= 0.5404 
poi = -0.09388 
3920 = -0.08935 
Pii =0.1549 
P02 = —0.35 1 5 
P30 = 0.005578 
P21 = -0.02021 

Pi2 = 0.05055. (1) 

Poo = 0.2273 
pio = -0.03987 
Poi = -0.02136 
P20 = -0.002942 

0.00655 
po2 = -0.001173 
3330 = 0.0007376 
P21 = -0.0005432 

Pi2 = 0.0001446. (2) 

Then, for M in Mj and S in ks per baryon, the following is an excellent fit to the R[S,M] function when 
R<2 Rj: 

LogR=poo + pwS + poiM + p2oS'^ + puSM + po2M'^ + p^oS^ + p2iS'^M + pi2SM'^ (3) 
R[S, M] » 10^°sR jij . (4) 



Likewise, for M > 2 Mj, 



